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Abstract 

The J/tp 7r — ► D D*, D D*, D* D* and D D cross sections as a function of 
y/s are evaluated in a QCD sum rule calculation. We study the Borel sum rule 
for the four point function involving pseudoscalar and vector meson currents, 
up to dimension four in the operator product expansion. We find that our 
results are smaller than the J/i/j ir — * charmed mesons cross sections obtained 
with models based on meson exchange, but are close to those obtained with 
quark exchange models. 
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Charmonium - hadron cross sections are of crucial importance in the context of quark- 
gluon plasma physics (![]. Small J /if) - hadron dissociation cross sections may favor an 
interpretation of the recent Pb + Pb data in terms of the production of a new phase of mat- 
ter. Part of these interactions happens in the early stages of the nucleus - nucleus collisions 
and therefore at high energies (y/s ~ 10 — 20 GeV) and one may try to apply perturbative 
QCD. However, even in this regime, nonperturbative effects may be important @. Inter- 
estingly, estimates using quite different methods give results clustering around the value of 
3 — 5 mb in this energy range. On the other hand, a significant part of the charmonium - 
hadron interactions occurs when other light particles have already been produced, forming 
a "fireball". Interactions inside this fireball happen at much lower energies (y/s < 5 GeV) 
and one has to apply nonperturbative methods. 

One possible nonperturbative reaction mechanism is meson exchange, which can be stud- 
ied by means of effective Lagrangians, constrained by flavor and chiral symmetries as well as 
by gauge invariance. This approach was first introduced in ref. [|3| and further developed by 
other groups Another reaction mechanism is quark interchange driven by Born-order 

matrix elements of the standard nonrelativistic quark model [pHTOj] . In this approach, once 



the masses and sizes of the mesons are fixed, there are no free parameters left. 

The results of the calculations for the charmonium-pion cross sections based on these 
two approaches can differ by two orders of magnitude in the relevant energy range. The 
situation clearly calls for different types of calculations that are constrained by other, in- 
dependent pieces of phenomenology. In this work we use the QCD sum rules (QCDSR) 
technique pT|JT5|| to study the J/if> — tt dissociation. The QCDSR technique allows one to 



compute hadronic quantities like masses, coupling constants and form factors in terms of 
quark and gluon properties and universal matrix elements which represent the properties of 
the QCD vacuum. In view of our relatively poor understanding of J /if) reactions in nuclear 
matter and considering the large discrepancies between different model estimates, we believe 
that our work adds to a better understanding of this important topic. 

We consider all four channels J /if) n — D D*, D D*, D D and D* D*. Let us start 
with the the four-point function for the process J /if) tt — > D D*: 



= i J d A x d 4 y d A z e~ ipi - x e ip3 - y e' 
x (Oim^jf (y)jJ(0)j D (z)}|0) , (1) 



with the currents given by j n = dij^u, j®* = u-y u c, = cry^c and j D = cij 5 d fTS] , where c, 
u and d are the charm, up and down quark fields respectively, and pi, p 2 , P3 and p 4 are the 
four-momenta of the mesons n, J /if), D* and D respectively, with p\ +p 2 = p 3 + p^. 

The phenomeno logical side of the correlation function, H^, is obtained by the consider- 
ation of J/if>, vr, D and D* state contribution to the matrix element in Eq. ([1]): 

U phen = m l F n mlfp m D *f D * m^f^p M af3 
^ m u + m d m c (p\ - ml) (p| - m 2 D ) 



x 12 Z^2 12 ZI2 H n. r. , [Z) 



9^a -P2^P2a/ml g v0 - p 3u p 3 p/ m 2 
p\ -m\ p\- m 2 D 

where h. r. means higher resonances and the hadronic amplitude for the process J /if) n 
D D* is given by 
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(3) 



We note that one has 1/pf pole in Eq. (0) in the limit of a vanishing pion mass. Fol- 
lowing p^-|15|, we can write a sum rule at p\ = and single out the leading terms in the 
operator product expansion (OPE) of Eq. ([1]) that match the 1 jp\ term. The perturbative 
diagram does not contribute with l/p\ and, up to dimension four, only the diagrams pro- 
portional to the quark condensate, shown in Fig. 1, contribute. After collecting the 1/p 2 
terms on the theoretical side and taking the limit pi^ —>■ in the residue of the pion pole, 
one obtains for the contribution of these two graphs 



r\<g<i> 



2m c {qq) pwipi^ + Pip, - 2p 3ll ) - px^v 



Pi 



(p 2 — m 2 c )(pl — m 2 c ) 



(4) 



Contracting the hadronic amplitude with the numerators of J/tp and D* propagators in 
Eq. (0) and comparing with Eq. (|J), the structure defining M.^ in Eq. (0) is easily identified. 
Therefore, defining 



= Add* (pi^Piu - Pi^u - 2p lu p 3/1 ) 



(5) 



we can write a sum rule for Add* i n any of the three structures appearing in Eq. (|^). To 
improve the matching between the phenomenological and theoretical sides we follow the 
usual procedure and make a single Borel transform, with all the external momenta (except 
Pi) taken to be equal: —p\ = —p\ = —pi = P 2 — > M 2 . The problem of doing a single Borel 
transformation is the fact that terms associated with the pole-continuum transitions are 
not suppressed [16| 



In ref. \\t 



it was explicitly shown that the pole-continuum transition 
has a different behavior as a function of the Borel mass as compared with the double pole 
contribution (triple pole contribution in our case) and continuum contribution: it grows 
with M 2 as compared with the contribution of the fundamental states. Therefore, the pole- 
continuum contribution can be taken into account through the introduction of a parameter 
Add* in the phenomenological side of the sum rule |13|,[14 ,[16|. Thus, neglecting m 2 in the 
denominator of Eq. 
get 



and doing a single Borel transform in — p\ = —p^ = —p\ = P 



wc 



A 



DD' 



+ A DD *M 



m 2 D , - ml 



-m? D /M 2 _ 



-m\IM 2 



m\-m 2 D 



-2m c (qq) 



-m 2 JM 2 



m c (m u + m d ) 



M 2 mlm 2 D mD*m i) F 7T f D fD*f,jj 



(6) 



where we have transferred to the theoretical side the couplings of the currents with the 
mesons, and have introduced, in the phenomenological side, the parameter Add* to account 
for possible nondiagonal transitions. 

At this point we should mention that the approximations we are using of exploiting 
the soft-pion limit and making the single Borel transform presents uncertainties. The main 
uncertainties are related to the continuum subtraction, the non-diagonal contributions, and 
the subtraction terms in the multiple dispersion relation. The approximation of the soft- 
pion limit can be ameliorated by going off the pion pole. In addition, further improvements 
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can be made by use of light-cone sum-rules [|L7|. These allow to use a light-cone pion 



distribution amplitude in substitution of condensates incorporating in this way additional 
QCD effects. Nevertheless, we believe that in this initial attempt our results should be useful 
as a comparison with what is obtained using model calculations. 

For consistency we use in our analysis the QCDSR expressions for the decay constants 
of the J/ip, D* and D mesons up to dimension four in lowest order of a s : 



3m 2 



du 



(u-m c ) (m 2 n _ u)/Mh 



mi 



m 



(qq)e (m D- m ^ /M M 



D 



fl* 



l7r 2 m 2 D , J,„- 



x e (™ 2 D * ~ S )/ M 



ds 



m. 



s \ s , 



(7) 



m D , 



■ (qq)e^- m ^ M M, 



f 2 = — 

H 47T 2 J Ami 



v, (r + 2m 2 ) Jr — Am 2 
dr e v '"^ 



3/2 



tJ^l-r)/M 2 M 



(8) 
(9) 



where represents the Borel mass in the two-point function. We have also omitted the 



numerically insignificant contribution of the gluon condensate. 

The parameter values used in all calculations are m u + m^ = 14 MeV 

m w = 140 MeV, m D = 1.87 GeV, m D * = 2.01 GeV, m^ = 3.097 GeV 

\3 n„\r3 



m r 



1.5 GeV, 
F n = V2f n = 

131.5 MeV, (qq) = — (0.23) 3 GeV 3 . We parametrize the continuum thresholds as um = 
(mM + A u ) 2 . The values of um are, in general, extracted from the two-point function sum 
rules for f D and f D * and in Eqs. (0), (|8|) and Using the Borel region 3 < M| f < 6 GeV 2 
for the D* and D mesons and 6 < < 10 GeV 2 for the J/ip, we found good stability for 
f D , f D * and with A u ~ 0.6 GeV. We obtained f D = 155 ±5 MeV, f D , = 195 ±5 MeV and 
= 225 ± 10 MeV, which are acceptable values for these decay constants [0. However, 
instead of using numerical values for these decay constants we use the sum rules in Eqs. (0), 
(0) and (P) directly when evaluating Ai. 



In Ref. 



it was found that relating the Borel parameters in the two-and three-point 
functions through M 2 = 2 M^, is a crucial ingredient for the incorporation of heavy quark 
symmetries, and leads to a considerable reduction of the sensitivity to input parameters, such 
as the continuum thresholds, and to radiative corrections. Therefore, we will use M 2 = 2 Mf^ 
to relate the Borel parameters and will work in the Borel range 8 < M 2 < 16 GeV 2 . We 
recall that this region corresponds to 4 < M 2 f < 8 GeV 2 , in which we have obtained good 
stability for the two-point sum rules of Eqs. (0), (|8|) and (^). This region also covers the 
range of the average values of the masses of the D, D* and J/ip mesons. 

In Fig. 2 we show, for A u = 0.6 GeV, the QCD sum rule results for Add* + A-dd*M 2 
as a function of M 2 (dots). We see that they follow a straight line in the Borel region 
8 < M 2 < 16 GeV 2 . The value of the amplitude A is obtained by the extrapolation of the 
line to M 2 = |T3|JT^JT6[| . Fitting the QCD sum rule results to a straight line we get 



A 



DD* 



17.71 GeV" 



(10) 
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As expected, in our approach A is just a number and all dependence of M.^ (Eq. @) on 
particle momenta is contained in the Dirac structure. This is a consequence of our low 
energy approximation. 

Next, we consider the process J /if) it — > D D (J/if> tt — > D* D*). In this case we have 
to change the current jf? (jo) in Eq. ([I]) to uijsc (cj a d). The phenomenological side is 
obtained as 

jjphen = m l F * m DfD 'day + P2aP2y 

M m u + m d m 2 (p| - m^) 

(Pi - ml) (pi - m^) (p| - m 2 D ) ' ' 

for J/^i 7r — > D D, where the hadronic amplitude is defined by M. = ■M- fJ ,ipi,P2,P3,P4) 
In the same way, for J /if) it — > D* D* we get 

nP he„ = mlF* m 2 D ,f 2 D , m^U 
^ va m u + m d {pi - ml) 

,2 



9tf + P2vP2p/m^ -g vS + p 3u p 3S /m 2 D , 



P\ ~ m \ P%- m l* 

x -9a<,+P^/mh +h r (12) 

with the corresponding hadronic amplitude defined by M. = -M ^uaipi, P2, P3, Pa) £2 e T e T- 
Similarly to the case J /if) it — > Z) £)*, in the OPE side the only diagrams, up to dimension 
four, contributing with l/p\ are the diagrams shown in Fig. 1. Therefore, taking the limit 
Pin — ► in the residue of the pion pole we get: 



n <«> = 2 (gg) ^aPxPzPl 
^ Pi (pi — m l)(pi — m c) 



(13) 



and 



> 



2{qq) 



pl(p 2 3 -m 2 c )(pj-m 2 

ml + P3-Pa) tanvpPx + , (14) 



where 



E^ ua — P\P\p\(—^vfi\^9ap, + ^^X-j9ap ~ ^aPX-yQ^u) 

+ e^pxipiPiPsa ~ P3P4PI0) + W/3a(PiPsPim 

+ P3PAP11* - PiPaPs» - PiPsP^) + ^fixi-PiPsPiv 

- P3PAPIV + P\P\P\v + PlPsPAfi) ■ (15) 

Comparing the phenomenological and OPE sides of the correlators we can identify the 
structure defining the hadronic amplitudes: 
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= A DD ^ a ^V\V2,v\ , M^ a = A D * D * E^ a . (16) 

It is important to notice that in writing Eq. fllED we have neglected the structure ea^ppf 
in M.^ ua . This is because, as can be seen from Eq. (|14D, this structure contains a term p^-Pi 
that can be rewritten in terms of p\ — m 2 and p\ — m 2 and, therefore, will contribute with a 
single pole which contains information about pole-excited states contributions. Since these 
contributions are considered in the phenomenological side as a parameter, we do not need 
to include them explicitly in the OPE side. 

We can write a sum rule for Add in the structure e^uPiPzPi, and a sum rule for Ad*d* 
in any of the structures appearing in Eq. ([15]). Thus, neglecting m 2 in the denominator of 
Eqs. (|TT| ) and ([12]), and doing a single Borel transform in — p\ = —p\ = —pi = P 2 , we get 

A MM + A MM M 2 2 m u + m d 

JM\. M ) — 



o 9 Jiw V 1 " / ^ 1V1 9 2 n / 



e 



-ml/M 2 



2 

where the subscript M stands for the D or D* mesons, with Cd = Cd* — 1 and 

rn D 

/M(M2) = 1_ - m ,_4 • (18) 

In Fig. 2 we also show, for A u = 0.6 GeV, the QCD sum rule results for Add + AddM 2 
(diamonds) and Ad*d* + Ad*d*M 2 (triangles) as a function of M 2 from where we see that, 
in the Borel region 8 < M 2 < 16 GeV 2 , they all follow a straight line. As explained before, 
the value of the amplitudes Add and Ad*d* are obtained by the extrapolation of the line to 
M 2 = 0. We get: 

Add ^ 12.25 GeV" 1 , Ad* d* ^ H-39 GeV" 3 . (19) 

Having the QCD sum rule results for the amplitude of the three processes J/ip tt — > 
D D*, D D, D* D*, given in Eqs. (|) and ([l|) with A given in Eqs. fllOD and (H) we can 
evaluate the differential cross section. 

Using our QCD sum rule result in Eqs. (|5]), (|T6|) , ( [TO]) and ( [19]) we show, in Fig. 3, the 



cross section for the J/ip 7T dissociation. It is important to keep in mind that, since our 
sum rule was derived in the limit pi —>■ 0, we can not extend our results to large values 
of t/s. Also, since the perturbative contribution is absent in our calculation, we were not 
able to properly disentangle the continuum contribution and our cross section may include 
contributions from higher states. Whereas they are certainly not important in the case of 
the pion, they may give some contribution to the heavy currents. Therefore our J/ip tt cross 
section may implicitly include (at least partially) the process ip' tt. For this reason, our 
numbers might be regarded as upper bounds. 

Our first conclusion is that our results show that, for values of \fs far from the J/ip tt — > 
D* D* threshold, (Jj/^^d*d* > °j/^->dd*+dd* > Vj/^^dd, hi agreement with the 
model calculations presented in |4j] but in disagreement with the results obtained with the 
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nonrelativistic quark model of ||, which show that the state D*D has a larger production 
cross section than D*D*. Furthermore, our curves indicate that the cross section grows 
monotonically with the c.m.s. energy but not as fast, near the thresholds, as it does in the 
calculations in Refs. 0-0]. Again, this behavior is in opposition to |J, where a peak just 
after the threshold followed by continuous decrease in the cross section was found. 

At higher energies, due to our low energy approximation, our approach gradually loses 
validity. In the fiducial region, close to threshold, 4.1 < y/s < 4.3 GeV, we find 2.5 < a < 
4.0 mb and these values are much smaller than those obtained with the effective Lagrangians 
without form factors in the hadronic vertices, but agree in order of magnitude with the quark 
model calculations of ||. 

Finally, we should mention that we have been studying the dissociation processes of the 
J/ijj in vacuum and the quantities relevant for QGP physics are in medium cross sections. In 
our approach the main effect introduced by the medium is the modification of the conden- 
sates, which is thought to be very mild. Our results depend only on the quark condensate 
and since it decreases with the nuclear density, we expect a further reduction in our cross 
section in a dense nuclear environment. 

In conclusion, we have used the QCD sum rule approach to evaluate the hadronic am- 
plitude of the J/tp tt dissociation. From the hadronic amplitude we have evaluated the 
J/ijj tt — > charmed mesons dissociation cross section, and have obtained 2.5 < a < 4.0 mb 
at 4.1 < t/s < 4.3 GeV. In view of the uncertainties discussed above these numbers should 
be taken as upper limits. 



It is interesting to remember that Bhanot and Peskin [2(J] have also used the OPE in the 
short distance limit to study the charmonium hadron cross section. This work was latter 
enlarged and updated by Kharzeev et al. [21] and also by Oh et al. [2^|. In these papers the 



crucial assumption was made that the charmonium is very small and resolves the partonic 
structure of the light hadron. In our approach we do not use this assumption, and we obtain 
larger values for the cross section. This seems to indicate that size effects are important, 
and that the J/ip cannot be considered as a nearly point like object. 

This work was supported by CNPq and FAPESP (contract numbers 98/06590-2, 
99/12987-5 and 00/04422-7). 
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FIGURES 




FIG. 1. Diagrams that contribute with l/p{, up to dimension four, in the OPE side of the 
amplitude vr + J/tp -> D + D* . 




FIG. 2. Amplitudes of the processes vr J/ip -> D D* + DD* (dots), D D (diamonds) and £>*£>* 
(triangles) as a function of the squared Borel mass M 2 . The solid, dotted and dot-dashed lines 
give the extrapolations to M 2 = (respectively). 
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FIG. 3. Total cross sections of the processes J/i/j it — > D D* + D D* (dashed line), D D (dotted 
line) and D* D* (dot-dashed line). The solid line gives the total J ftp it dissociation cross section. 
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